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ABSTRACT 

Well-calibrated spectropolarimetry studies at resolutions of R >10,000 with signal-to-noise ratios (SNRs) better than 0.01% across 
individual line profiles, are becoming common with larger aperture telescopes. Spectropolarimetric studies require high SNR obser¬ 
vations and are often limited by instrument systematic errors. As an example, fiber-fed spectropolarimeters combined with advanced 
line-combination algorithms can reach statistical error limits of 0.001% in measurements of spectral line profiles referenced to the 
continuum. Calibration of such observations is often required both for cross-talk and for continuum polarization. This is not straight¬ 
forward since telescope cross-talk errors are rarely less than ~1%. In solar instruments like the Daniel K. Inouye Solar Telescope 
(DKIST), much more stringent calibration is required and the telescope optical design contains substantial intrinsic polarization ar¬ 
tifacts. This paper describes some generally useful techniques we have applied to the HiVIS spectropolarimeter at the 3.7m AEOS 
telescope on Haleakala. HiVIS now yields accurate polarized spectral line profiles that are shot-noise limited to 0.01% SNR levels 
at our full spectral resolution of 10,000 at spectral sampling of ~ 100,000. We show line profiles with absolute spectropolarimetric 
calibration for cross-talk and continuum polarization in a system with polarization cross-talk levels of essentially 100%. In these data 
the continuum polarization can be recovered to one percent accuracy because of synchronized charge-shuffling model now working 
with our CCD detector. These techniques can be applied to other spectropolarimeters on other telescopes for both night and day-time 
applications such as DKIST, TMT and ELT which have folded non-axially symmetric foci. 

Key words. Instrumentation: polarimeters - Instrumentation: detectors - Techniques: polarimetric - Techniques - spectroscopic - 
Methods: observational 


1. Introduction 


Polarimetry and spectropolarimetry generally require high 
signal-to-noise observations from intrinsically imperfect tele¬ 
scopes and instruments. The fundamentally differential nature 
of polarimetric measurements is an advantage for minimizing 
systematic instrumental noise. Nevertheless, achieving photon- 
noise limited accuracy is often dependent on mitigating subtle 
instrument systematics. These errors fall into two categories: I) 
calibration uncertainty, and II) measurement instability. We de¬ 
scribe here a strategy for minimizing both of these uncertainties 
using new technologies and algorithms in the context of stellar 
spectropolarimetry. 


Measured optical polarization is sensitive to the local geom¬ 
etry of the radiating or scattering source and allows otherwise 
optically unresolved features of the source to be inferred using 
forward modeling approaches. We have been investigating po¬ 
larimetric signatures in individual spectral lines of young stars 
that show Stokes QU features with amplitudes of about 0.01% 
(iHarrington et al.ll20T^ : lHarrington & Kuhnll2009bll2008l) . Such 
tiny polarized spectral features can contain basic information 


about, for example, the geom etry of the circumstellar environ- 
ment (iKiihn et al .11201 lLl2007h . 

This problem of measuring small polarized features in astro¬ 
nomical objects is not at all unique to stellar spectropolarimetry 
and has a long heritage beginning in solar magnetic field mea¬ 
surements. Most measurement techniques that deliver such high 
precision usually follow a multi-step process. Some type of mod¬ 
ulation occurs where polarization signals are optically converted 
to intensity variations. To de-modulate these variations, several 
independent exposures are somehow differenced to remove some 
systematic errors. These differential measurements then have to 
be calibrated to ultimately derive an absolute polarization of the 
input signal. Several calibration techniques are usually required 
to fully recover the input polarization state of the source from 
the polarizing influence of the modulators, telescope and detec¬ 
tors. One important calibration need is to correct the scrambling 
of the input polarization state from the observed output polariza¬ 
tion state (called cross-talk), that is caused by the telescope and 
detectors. At a lower level telescopes also often create or destroy 
the degree of polarization of the input light. This induced po¬ 
larization and depolarization, must also be calibrated. Different 
modulation strategies, beam differencing techniques and techni- 
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cal considerations drive instruments to very differ(^nt approaches 
depending on the science goals and practical limitations. For 
some specific i nformation on typical n iodulatior schemes an d 
efficiencies , see iKeller & SnikI (l2009l) : IXomczvkTt al.l (l201Qh : 
IWiin et ^ (1201 ll) , astronomic al measuremen t techniques and 
specific stellar applications, see lSnik & Kellerl T^. 

Calibrating night-time spectropolarimetry gen orally requires 
data from both known polarized and unpolarized sources. Other 
techniques make assumptions about for example the polarized 
spectral symmetry of the source; or partial knowliedge about the 
polarizing properties of the telescope and instrum mt. For exam¬ 
ple, solar observations can be calibrated to a high level of ac¬ 
curacy using observations of s ources with a fund mental wave¬ 
length polarization symmetry (iKuhn et al.lll994 ). In all cases 
the full Stokes spectra must be measured in order to achieve 
the most precise absolute polarization calibration. Under certain 
limitations and assumptions, the number of meamrements can 
be reduced for efficiency reasons. For example, instruments with 
relatively low cross talk can neglect to record sor le parts of the 
Stokes vector if they are deemed unnecessary or i lefficient. 

Minimizing type II ("stability") polarization uncertainty re¬ 
lies mostly on fixing instrument, telescope, and atmosphere 
variability during the modulation process wher^ differencing 
intensity measurements is necessary to derive tie full-Stokes 
polarized spectra. Techniques like dual beam analyzers, fast- 
modulation, beam swapping modulation, or other instrument 
stabilizing schemes decrease this noise in sp(^ctropolarime 
ter-rtelescope measurements. Where many photons are needed, 
obtaining efficient calibrations is especially impo lant. The ma¬ 
jor limitations are often observing duty-cycle 
drifts, waste of precious night time to calibrate, 
to calibrate at all. We demonstrate here a new instilument using a 
combination of techniques that helps minimize s(^veral of these 
error sources at a highly non-optimal polarimetric telescope fo 
cus. With these techniques we can deliver both sta 
differential accuracy with reasonable absolute po^ 
bration without rebuilding the telescope to be optimized for po- 
larimetry. 


instrumental 
3r an inability 
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1.1. Spectropolarimetric instruments, errors & Suppression 
techniques 


d^ 


Large aperture telescopes are tasked for various p^ 
quirements. For example, some are optimized to 
continuum polarization, low cross-talk, or high s 
tion and sensitivity. On short timescales telescof 
ing errors and atmospheric seeing can change tl 
a beam on the telescope and spectrograph optics 
ulate the detected intensity to create sp urious 


pci; 


long timescales optical coatings oxidize (Ivan Hart cn et al.ll20Q9l) 


making calibrations unstable in time. As non-e ijuatorial tele¬ 
scopes track a target, mirrors rotate with respect to the field to 
change the properties of the optical system. Scattered light and 
instrument instabilities in a spectrograph can cause drifts. The 
electronic properties of the detector pixels, readout amplifiers 
and signal processing algorithms imprint fiat-fieldi ng errors, bias 
drifts and other non-linearities into the data. Complex optical 
paths such as coude foci involve fold mirrors at oblique angles 
that induce strong polarization effects. 

Many of these errors can be overcome by usi ng calibration 
techniques or by designing special-purpose instru nents. Formal 
polarimetric error budgets applied over wide ran ^es of science 
cases are typically not imp lemented but formalis r is to do so are 
becoming more common (iKeller & SnikI l2009l : Ide Juan Qve- 
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record two orthogonally analyzed beams in a single 
and double the instrument efficiency. Most instruments 
-exchange or other modulation techniques to difference 
t exposures and remove many effects caused by detec- 
nics, fiat-fielding and optics behind the analyzer, 
ost all dual-beam night-time spectropolarimeters, sev- 
lendent exposures are required to make a full-Stokes 
Jnent. The timescale for these exposures (and polari- 
odulation) is often many minutes. For instance, there 
existing and planned fiber-fed instruments like ES- 
HARPSpol, BOES, PEPSI that are designed to deliver 
ijesolutions above 50,000 an d large one-shot wav elength 
by using cross-dispersion (ISemel et~^ll993l: Donati 
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ISnik et 
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et al. l20C 8h. These instruments are not designed to measure con- 
p Dlarization but fairly reasonable m easurement precisio n 
e under a range of circumstancesFI (iPerevra et al.l2015l) . 
these instruments do not have spatial resolution due 
er-feed. In addition, absolute calibration can be lim- 
oss-talk systematic errors in up-stream optics even at 
n focus. For ESPaDOnS, several years of observations 
time-dependent cross-talk errors from its atmospheric 
compensator and collimating lens. These effects lim- 
ute calibrations to the >3% level in the initial few years 
Ijnproved optics and testing reduced the cross-talk to be- 
% design thr esholds, but it is still present as the major 

_ n limitation ([Barrick et af]|201Ql : iBarrick & Benedict 

l201Ql : IPe ’evra et al.ll2015l) 
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1201 2l : iDe Juan Qvelar et al.ll2012l) . Dual-beam po- 


1999l:lPiskunov et al.l201ll:lManset & ponatill2003l: 

1l201ll:lKim et al.ll2007 : ISnik et al.ll2008l: Strassmeier 


Othei • night-time slit-fed instruments have different perfor- 
ipabilities and correspoi idingly different stab ility is- 
as FQRS on the VLT (iBagnulo et al.ll2012l: Seifert 


d: IBagnulo et al.l 120021 l2009l: iPatat & Romaniellol 

assegrain instruments such as LRISp on Keck, ISIS 
HT are examples of variable-gravity spectropolarime- 
h are designed to deliver relatively stable continuum 
on, imaging capabilities and more modest resolutions 
can be subject to mechanical flexure or environment 
unless designed and stabilized properly ( Goodrich & 


Goodrich et al.lll995l: lLeonell2007l: Biork man & 
Lomax et al.ll2012l: iNordsieck et al.ll200ll: Nord- 


I Iarris ll996l : lEversberg et al.lll998l) . 
fen true that cross-talk errors are strong and time- 
in instruments where polarimetry was not a driv- 
itement. Many large night-time telescope projects (such 
AT, GMT or ELT) do not allow access to Cassegrain 
ian foci for traditional polarimetric instruments with 
talk. Some of the most accurate single line spec- 
netric measur ements ever obta ined have been realized 
is telescopes (iLin et al.ll^04h . The great advantages 
ic range and scattered light stability of such tele- 
1 itweighs the fractional percent polarization induced by 
metric primary optics. Eor instance, the Daniel K In- 
r Telescope (DKIST) will make good use of several 
lution spectropolarimeters operating at the 10“"^ pre- 
els. This system relies on several precise calibration 
s that impart stringent requirements on the telescope 
Mueller ikatrix (iKeil et al.l201 ll:lKelled2003l:lNelson et al.l20ia: 


et al .11^121: ISocas-Navarro et alJ ^05h. As projects 


grow m fecope and complexity, real error analysis and tradeoff 
studies are becoming essential as exemplified by the VLT and 


^ http://ww w.cfht.hawaii.edu/Instruments/Spectroscopv/Espadons/ContiPolar/ 
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E-ELT codes under develop ment (Ide Juan Qvelar et~^l2012l : 
iDe Juan Ovelar et'n]l2012l) . Even axially symmetric telescope 
beams require calibration of instrumental polarizat ion from the 
instrument optics and non-uniformities in coatings dHough et aLl 


I 2 OO 6 I : IB ailev et~n]l20 1 OL l2008l : IWiktorowicz & Matthewdl2008h 

In the solar community, several instruments have included 
very fast polarimetric modulation to overcome many of these 
common instrument issues to achieve 0.001% differential pre¬ 
cision polarimetry. Modulators include liquid crystals (LCs), 
piezo-elastic modulat ors (PEMs), rap i dly rotating compo nents 
and other te chniques (lGandorfedll999l:lElmore et^ll992l: Sku- 
manich et al. ll997l : lGandorfer et al.ll2004lLin et al.ll2004l) . Mod- 
ulation can be chromatic ally balanced, tuned or optimip d for 
various o bserving cases (|Povellll995l: iGisler et al.ll2003l: Tom- 
czyk e t al. l2010l:lde Wiin et al .1120101: 1 Wiin et al.ll201 lUSnik et al.l 
l2009l : iLopez Ariste & Semel l201ll : iNagaraiu et al.l 2007h . At 


modulation rates faster than IkHz most instrument effects and 
even atmospheric seeing fluctuations can be suppressed ( Keller 
et al. 11994 IStenflol 120071: iHough et al.ll2006l : [Hanaoai2004 
iRodenhuis et al.l1^l2l : IXu et al .I l2006l) Night-time modulation 
frequencies are usually much slower and fixed at the exposure 
time so that detector read noise does not dominate. 

Other strategies to calibrate, stabilize and optimize spec- 
tropolarim eters with varying degrees of cros s -talk have been de ¬ 
scribed (cf.lElmore et al. (|2Q10l): Snikl(l2006l):lTinbergen (2007); 


Socas-Navarro et all (l201 lh:ldel Toro Iniesta & Collados (2000) 


Tinberpij (l2007h : IWitzel et^ (I2011 ) : iBagnulo et al. (I2009L 

2012h : lGiro eTall ((2003l) : IWitzel et all (l201lh : lMaundl (l2008h ). 


Some use spectral modu lation or other te chniques to accomplish 
instrument stabilization (ISnik et al.ll2009l) . Other polarimeter de¬ 
signs, such as the VET X-Shooter have controlle d chromatic 
prope rties and use special techniques to calibrate (ISnik et al.l 

I2012h . 

In night-time spectropolarimeters, magnetic field studies of 
stars are pushing precision limits by co-adding many lines into 
a single pseudo spectral line that yields effective signal-to-noise 
ratios above 100,000. Large surveys such as MIMES 0 and oth¬ 
ers use ESPaDOnS and NARVAL fiber-fed spectropolarimeters 
at sp ectral resolutions greater than 60.000 are exam ples (Wade 
et al. I 2 OIOI : iGrunhut et al.ll2012l : IWade et~.ll2012l) . These pro¬ 
grams combine many individual spectral lines over many spec- 
tral orders to derive a single representative pseudo line profile 
(iDonati et all 1 199^: iKochukhov et al1l2010l: j^ennhauser et all 


120091: IKochukhov & Siidnikll2ni 3l:lPetit et al .112011h . 

In other night-time science applications, single-line studies 
at more modest resolutions of 3,000 to 10,000 are sufficient 
but require similarly high signal-to-noise ratios across individ¬ 
ual lines. Eor instance, Herbig Ae/Be stars and other emission¬ 
line stars show polarization signatures at the 0.1% level and be¬ 
low. However, individual lines have varying profiles and line- 
combination techniques don’t apply and sufficient SNRs are 
require d to see individua l line profiles (iHarrington & Kuh^ 
l2nn9hini2(Ml2nn7Ll2ni 1h . Herbig Ae/Be and T-Tauri polarized 
line profiles can be quite large and variable when obse rv ed wit h 
effective resolutions of even 3,000 (IVink et al.ll2005allRl2002l) . 
Absolute continuum polarization is a particularly important mea¬ 
surement for Be star and some other emission-line studies be¬ 
cause interpreting the polarizing mechanisms depends on know¬ 
ing pol arization ciirection t o distinguish models fo r the source 
region (lBiorkmanll2006allE: iBiorkman et al.ll2000allE: Brown & 


McLean ll977l:ICarciofi et al.H2009l:lHarries & Howard 


- ^ 

1 19961 : llgnacel200Ql : lignace et ^1 19991 : lOudmaijer & Dre^ 


It] 


1997 


1999 


lOudmaiier et al.ll2005l : lOuirrenbach et al.l 1997 ). Wolf-Ray e t 
stars can be similarly probed (iHarries et al.ll200(]l:IVinkll201(]h . 
Its important to note that there are a large number of interesting 
stellar problems where spectropolarimetric resolutions as low as 
3000 are sufficient to learn about the source region. 

We upgraded AEOS HiVIS spectropolarimeter to minimize 
several sources of instrumental instabilities to achieve high 
SNRs and reasonable calibration accuracies by using liquid crys¬ 
tal modulation and charge-shuffling detectors. HiVIS has a spec¬ 
tral sampling of 4.5pm to 5.5pm per pixel corresponding to a 
1-pixel sampling of R~ 150,000. Eor this campaign we use a 
wide slit delivering a spectral resolution of 12,000 which is over¬ 
sampled by a factor of about 9 when charge shuffling. SNRs 
across the Hq, line profile of 0.01% are now achievable in this 
mode on fairly bright stars in a single night while delivering ac¬ 
curate continuum polarization with a stable polarization cross¬ 
talk calibration. We have implemented faster modulation tech¬ 
niques combined with electronic charge-shuffling to achieve this 
instrument stability. We have also improved calibration tech¬ 
niques using the daytime sky, demodulation techniques and other 
polarimetric optics in order to achieve reasonable continuum po¬ 
larization stability and absolute polarization calibration. 

We show here that night-time long-slit spectropolarimetric 
observations can be calibrated for cross-talk to an absolute po¬ 
larization reference frame to roughly ~ 1 % levels in the presence 
of 100% cross-talk. The polarization calibration changes with 
wavelength are generally smooth functions that can be measured 
at much lower spectral resolution. When calibrating spectral line 
profile measurements with 0.01% signal-to-noise limits, these 
absolute calibrations can be treated simply as constants. Other 
technologies we use allow us to stabilize the stellar continuum 
polarization to a 0.1% calibration level. These combined tech¬ 
niques allow us to deliver both high differential precision across 
line profiles while maintaining good absolute polarization cal¬ 
ibration even in the presence of several difficult instrumental 
challenges. 

1.2. HiVIS and the charge-shuffling spectropolarimeter 

We have been upgrading the High Resolution Visible and Infra¬ 
red Spectrograph (HiVIS) on the 3.67m Advanced Electro- 
Optical Syst em (ALPS) telescope on Haleakala. Maui (Har¬ 
rington et al.l20 id: IHarrington & Kuhnll2008l:lHarrington et al 


2006 : lHarringtonl2008l : IHarrington et al.ll201irrThornton et al 


http://www4.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/MiMeS/ 


2003I) . HiVIS is a long-slit spectrograph built for spectral reso¬ 
lutions of 10,000 to 50,000 in the AEOS coude room. HiVIS is 
a cross-dispersed echelle spectrograph designed to deliver com¬ 
plete wavelength coverage with high order overlap across 15-20 
spectral orders on a 4k by 4k focal plane mosaic. An optical lay¬ 
out of the spectrograph ^d polarimeter can be found in H arring - 
ton et al. (l2006l) . lThornton et al.l (l2003h and lHarringtonl (l2008h . 
Eigure \T\ shows a schematic layout of the telescope and spec¬ 
tropolarimeter. There are 8 optical elements before the beam is 
delivered to the coude room. The primary and secondary mirror 
create a converging f/200 beam after the secondary mirror reflec¬ 
tion. This beam is delivered to the coude room by 5 fold mirrors 
at 45° reflection angles. There is also a selectable window of 
either BK7 or Infrasil between the last two fold mirrors. After 
entering the coude room, there are 8 more reflections before the 
slit and the polarimeter in this HiVIS charge-shuffling configu¬ 
ration. These optics accomplish collimation, tip/tilt steering at 
a pupil image and basic optical packaging. With this complex 
optical layout, our instrument causes linear-to-circular cross talk 
of 100% at some telescope pointings and wavelengths. At the 
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Fig. 1. This Figure shows a schematic layout for the AEOS and HiVIS optical pathway when in the charge-shuffling configuration. The telescope 
optics are shown in blue. There are 5 fiat fold mirrors that bring the beam in to the coude room. The elevation and azimuth axes are labeled to show 
the optical changes during telescope tracking. The HiVIS optical relay which feeds both visible and infrared spectrographs and accomplishes tip/tilt 
correction is shown in green. The spectrograph optics are shown in red. All optical surfaces are shown. Polarization and calibration components 
are shown in black. The HiVIS modulator and analyzer are behind the spectrograph slit. All optics in blue and green contribute to telescope 
polarization cross-talk. The Full-Stokes Injection Unit in front of the HiVIS spectrograph slit is for polarization calibration of the modulators and 
spectrograph. Additional calibration optics shown in black between the AEOS telescope (blue) and relay optics (green) are mounted at the entrance 
of the coude room. These provide fiat fielding, wavelength and polarization calibration light sources. See text for details. 


entrance to the coude room we have a calibration unit with arc, 
and flat field lamps.. There is also a polarization calibration unit 
consisting of an achromatic quarter wave plate and a wire grid 
polarizer under computer controlled rotation. The AEOS tele¬ 
scope has the largest telescope polarization we have experienced 
and is a challenging system to calibrate. 

We describe here the charge shuffling and liquid-crystal 
modulator improvements and the calibration algorithms we use 
for obtaining accurate full-stokes spectropolarimetry of stellar 
and circumstellar, visible and near-IR sources. These techniques 
yield continuum polarization stable to <0.1% accuracy across all 
spectral orders. 

1.3. HiVIS Polarimeter Upgrades 

In order to accomplish liquid crystal modulation and syn¬ 
chronous charge shuffling two tasks had to be accomplished: 

1. Enable clocking of photoelectrons under software control 
during an exposure. 

2. Synchronize charge clocking with liquid crystal modulation 
state control. 

This allows us to record the 4 spectra (2 modulation states in 
a dual-beam optical path) in one exposure. The detector control 
electronics were upgraded using the STARGRASP hardware and 
software package. This package was developed for orthogonal 


transfer aPDlications (iTonrv et al.lll997l:lBurke et al.ll2007l: Qn- 
aka et al. l2008l) . This new controller package will move accumu- 
lating photoelectrons an arbitrary number of pixels and arbitrary 
number of times using software commands during an exposure 
and before detector readout. The charge packet accumulating in 
each pixel of a CCD column can be clocked both towards and 
away from the shift register with simple command-line tools. 
The typical row transfer time depends on several settings but is 
less than lOy/s for our devices. Since the focal plane of cross- 
dispersed echehe instruments naturally have separated spatial 
and spectral directions to the focal plane illumination pattern, 
this clocking is set to be along the spatial direction of the slit 
image which is also the readout direction. 

Normal CCD’s typically do not allow charge transfer in 
two orthogonal directions but our controller firmware can clock 
charge away and towards the shift register. By synchronizing 
these shifts with changes in modulator states, adjacent pixels 
can be uses as storage buffers without adding detector readout 
noise. If the order separation is large, the charge can be electron¬ 
ically smeared in the spatial direction to effectively increase the 
full-well capacity during the exposure. This allows much more 
efficient use of the pixels in the CCDs. 

Meadowlark nematic liquid crystal modulators (LCs) were 
installed and a basic polmzation calibration analysis was re¬ 
ported in iHarrington et al.l (|201Q|) . The HiVIS control software 
was updated to include synchronous command of the modulators 
and the charge shuffling on the focal plane mosaic. The STAR- 
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Fig. 2. This Figure shows examples of data delivered by the different modes available to Hi VIS users. Each row shows a different observing 
mode. The left column shows a small portion of a single exposure with the wavelength direction as generally horizontal and the spatial direction 
as vertical. The normal dual-beam polarimetric mode is shown in the top row. Two orthogonally polarized orders are imaged on the CCD. The 
corresponding spatial profiles are computed as a spectral average in the spatial direction across an order. The spatial profile is shown on the top 
right. Two beams are seen with a FHWM of roughly 10 pixels. The Savart plate analyzer separates the two orthogonally polarized beams by roughly 
54 pixels aligned with the spatial direction of the spectral orders. The dekker is set to a slit length less than half this displacement in the spatial 
direction. The liquid-crystal charge-shuffled mode is shown in the middle row. The two distinct liquid crystal modulation settings are inter-woven 
in 4 beams recorded in a single exposure. We use a charge-shuffling distance of 27 spatial pixels, half of the Savart plate beam displacement. Two 
beams are stored for one modulation state in non-illuminated buffer pixels while the photoelectrons are accumulated for the second liquid crystal 
modulation state. Four beams are recorded using two illuminated spectral regions and adjacent non-illuminated buffer regions. The smear-mode is 
shown in the bottom row. The accumulating charge packet from both orthogonally polarized beams is slowly shuffled to near-by non-illuminated 
rows as the exposure progresses. This effectively increases the full-well capacity of the image and increases the saturation limit of the detector. 


GRASP package uses ethernet communications over a local gi¬ 
gabit network to synchronize external devices with the charge 
shifting. The timing jitter is less than 10 milliseconds for our 
ethernet commands. The LCs have a nominal switching time of 


40 milliseconds. Thus timing jitter between the charge shuffling 
and the liquid crystal switching will be small and will typically 
be averaged over many switching cycles. 
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The last optical requirement was to substantially reduce the 
slit height with a dekker so the 4 recorded beams were optically 
separated by our Savart plate analyzer. A dekker width of 30 
spatial pixels or less must be used with our Savart plate to en¬ 
sure optical order separation. The Savart plate displacement is 
between 53 and 55 pixels thus we charge shuffle by 27 pixels in 
the default operating modes. 

1.4. HiVIS Spectrograph & New Observing Modes 

With the new 4k by 4k mosaic focal plane, we can continuously 
cover a wavelength range >220nm typically in the 550nm to 
950nm range with good overlap between spectral orders. The 
order overlap is usually between 10% and 40% depending on 
the order and cross-disperser settings. The optical design deliv¬ 
ers a two-pixel per 0.35“ monochromatic slit image for the high¬ 
est spectral resolution (R -50,000) settings. This requires spec¬ 
tral sampling in the range of 4pm to 6pm per spectral pixel. For 
charge shuffling at high signal-to-noise ratios, we adapted the 
dekkar to use the 1.5“ slit for this campaign with a spectral reso¬ 
lution of R =12,000. We strongly over-sample in the spectral di¬ 
rection with sampling of roughly 4.3pm at Hq,. This allows us to 
efflciently observe bright targets to high SNRs without saturation 
while minimizing several kinds of instrument instability errors. 
Note ESPaDOnS spectra have quite similar spectral sampling to 
HiVIS. ESPaDOnS samples at 2 points per monochromatic slit 
image for R=68,000 and we sample at 2 points per R=50,000 
monochromatic slit image. Since the HiVIS slits are selectable 
with a filter wheel and our typical science cases do not require 
high spectral resolution, we used a wide slit for high throughput 
at R= 12,000 while maintaining the same spectral sampling. 

The modulator mechanical package is flexible and can be 
configured for liquid crystals or rotating achromatic retarders 
in a few minutes. Additionally, the HiVIS control software has 
three main observing modes. The first is standard dual-beam po- 
larimetry which only records one modulation state in two beams. 
The second mode is the liquid crystal modulated charge shuffling 
technique where four independent spectra are recorded, two for 
each modulation state. The third is called smear-mode where the 
accumulating charge packet is slowly shuffled back and forth 
during an exposure to more evenly distribute the charge in all 
available spatial pixels. This smear mode is used in bright star 
campaigns to effectively increase the saturation limit and reduce 
readout overheads by more fully using the available full well 
capacity of all available spatial pixels. Figure [2] shows small re¬ 
gions of an image for each observing mode on the left and the 
corresponding spatial profile on the right. 

There are several advantages gained by using these observ¬ 
ing modes for high signal-to-noise observations. Typically bright 
targets saturate the detector and require exposure times substan¬ 
tially shorter than the readout time. By implementing charge 
shuffling we double the time spent integrating for a fixed readout 
time and substantially increase observing efficiency. 

In the case of long-slit instruments, the delivered spatial pro¬ 
file typically takes on the shape of the atmospheric seeing com¬ 
bined with any beam jitter and tracking errors. This profile is nor¬ 
mally strongly peaked and Gaussian-like in shape where the core 
saturates leaving many underexposed spatial pixels. By charge 
shuffling, this accumulating charge can be more evenly dis¬ 
tributed across the focal plane. Beam location changes induced 
by the guiding and tracking system errors also cause changes 
in calibration and systematic uncertainties in several parameters. 
By doing a full dual-beam modulation sequence within a single 
exposure, many of these effects are mitigated. The liquid crys¬ 


tals are easily tuned via software commands to have different 
modulation schemes. 

The calibration of the liquid crystals is easily accomplished 
with our polarization calib r ation u nit shown in Figure [T] and out¬ 
lined in iHarrington et al.l (l201Qh . Additionally, we installed a 
new dichroic on our slit-viewing guider camera. This allowed 
faster and higher quality tip/tilt correction using a pupil steering 
mirror at IHz rates for targets brighter than V=8 under nominal 
conditions. 
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Fig. 3. This Figure shows the flat field lamp intensity spectrum (blue) 
and polarization errors (black and red) for a single spectral order. We 
show the data at our nominal 12,000 spectral resolution and R= 100,000 
spectral sampling. The right-hand y axis shows the intensity spectrum 
in detected counts per exposure. The intensity spectrum in blue shows 
a strong blaze function intensity variation across the order along with 
the typical CCD mosaic gaps and pedestal offsets. The two CCDs in 
the mosaic each have two amplifiers with slightly varying gains and bi¬ 
ases. The intensity was left uncalibrated to illustrate these CCD effects. 
Spectral features from the lamp and atmosphere are seen in the middle 
of the spectrum. The left-hand y axis shows the polarization spectrum 
multiplied by 10,000. Amplifier and CCD gaps are not visible in the 
polarization spectrum. The computed polarization is dominated by shot 
noise below the 0.01% level at full spectral sampling. The RMS of the 
polarization spectrum is substantially decreased by averaging adjacent 
spectral pixels, as seen in the red curve. This averaging shows that pho¬ 
ton statistics are the dominant noise source. See text for details. 


1.5. Liquid Crystal Modulation Rates 

To determine the optimum liquid crystal modulation rate, the de¬ 
polarizing effect caused by the nominal 5% to 95% liquid crystal 
retardance switching time of 40 milliseconds speed was tested. 
Figure ID shows a few select modulation matrix elements derived 
as functions of the switching time when using our polarization 
calibration unit at the HiVIS slit. As expected. Figure [4] shows a 
general reduction in modulation matrix element amplitudes at 
switching times faster than 500ms (or a 2Hz modulation fre¬ 
quency). As the liquid crystals are driven faster, they spend a 
proportionally longer time switching than resting, causing a de¬ 
polarization of the input signal. For our typical use case, mod¬ 
ulation of 2Hz or slower is selected. Under these conditions the 
number of cycles is typically low (e.g. 120 shuffles for a 60 sec¬ 
ond exposure). 
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Fig. 4. The modulation matrix elements as the LCVR and charge¬ 
shuffling system is driven at different shuffling speeds. The Meadowlark 
LCVRs have a nominal switching time of 20-50ms depending on the 
temperature, wavelength and retardance change commanded. As the 
modulation speed is increased the liquid crystals spend a larger frac¬ 
tion of the integration time in intermediate states while switching retar¬ 
dance, thus reducing the efflciency of the polarization modulation. We 
nominally use a switching time of 1 second or slower to achieve high 
polarization modulation efflciency with Hi VIS. 


We have extended calibrations done with LoVIS to the new 
Hi VIS charge-shuffled liquid crystal modulated configuration. 
As a standard non-precision tuning, we set our liquid crystal 
voltages to those derived in the laboratory for nominal Hq, mod¬ 
ulation using a Stokes definition scheme where each exposure 
measures a single Stokes parameter QUV. With the polarization 
calibration unit, we inject six pure + quv Stokes parameters we 
can derive two redundant sets of demodulation matrices. With 
this information, we also can extract the chromatic effects of the 
calibration quarter wave plate as well as the stability of the sys¬ 
tem. 

The derived liquid crystal modulation matrices are quite sta¬ 
ble over the course of a night. Each night we observe, we collect 
two complete modulation matrices derived using the -f and - in¬ 
put states. Experiments done many times show that independent 
sets of calibration observations can be used to demodulate other 
sets of observations to within the shot noise of <0.001 in any 
individual Stokes parameter. When a demodulation is done us¬ 
ing the same input state, say -f Stokes inputs demodulating mea¬ 
surements of -f Stokes inputs, the residuals are always shot-noise 
limited. When demodulation is done using opposite input states, 
the dominant effect is the chromatism in the calibration quarter 
wave retarder. This effect is removed when using the daytime 
sky calibrations to calibrate the cross-talk in the demodulated 
spectra. 

1.6. Polarimeter Summary 

With these optical improvements, we performed several tests to 
demonstrate the performance of this new hardware. In summary, 
the liquid crystals achieve >95% modulation efflciency at speeds 
of IHz or slower and the calibrations are stable over several-day 
timescales. The scattered light, bias stability and calibration data 
all give consistent results. In a series of 400 exposures, we were 
able to demonstrate the repeatability of the polarization calibra¬ 
tions down to the 0.01% SNR limit at full spectral resolution 


across the spectrum in the presence of molecular band features. 
Eigure[3]shows polarization calibrations where a set of 200 im¬ 
ages was used to calibrate another 200 independent Stokes q ex¬ 
posures. The results are stable with a signal-to-noise ratio (SNR) 
greater than 10,000 at full spectral sampling and resolution. Eig- 
ureOshows the intensity spectrum across a single spectral order 
along with the corresponding residual error between the two in¬ 
dependent 200-image calibration sets. There is a strong blaze 
function and some spectral absorption features from the lamp 
and atmosphere seen. Even with these intensity variations, spec¬ 
tral smoothing does not reveal underlying systematic effects that 
limits the instrument. All spectral orders show only statistical 
noise when spectrally binned by factors up to 128. The corre¬ 
sponding calibrations obey a ^fN behavior down to 0.001% lev¬ 
els. As an example in Eigure[3l the red line shows a simple box¬ 
car smoothing and a corresponding reduction in statistical noise 
of the spectra. Liquid crystal instabilities, calibration inaccura¬ 
cies or other detector-based effects do not appear. We conclude 
that our coude instrument does not have substantial systematic 
errors that limit the calibration accuracy. Our instrumental arti¬ 
facts are below 0.01% levels at full spectral sampling across the 
full spectral range for internal calibration sources. 

2. Calibrations: Steiiar Standards 

Observations of unpolarized and polarized standard stars are 
used to derive instrument calibration stability and precision lim¬ 
its. With the ability to record 4 modulated beams in a single ex¬ 
posure with IHz modulation rates, we can reduce the continuum 
polarization errors caused by instrument instabilities and atmo¬ 
spheric changes on timescales of seconds to minutes. 
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Fig. 5. This Figure shows the computed polarization spectra for re¬ 
peated measurements of the unpolarized standard star HR2943. Slit- 
guiding and charge-shuffling was disabled to demonstrate the typical 
effects introduced by atmospheric transmission variation (clouds), tele¬ 
scope jitter, guiding errors and other instrumental errors on the con¬ 
tinuum polarization. The intensity spectrum in the bottom row shows 
all 19 extracted spectral orders and the associated blaze functions. Each 
quv continuum polarization measurement shows substantial variation of 
over 2%. There are two components to the quv spectra: the spectrograph 
calibration function and the telescope induced polarization. See text for 
details. 

HiVIS observations of unpolarized standard stars without 
guiding or charge shuffling have shown the derived continuum 
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polarization to have high variance of up to 5%. Several un- 
polarized standyd star lists a re ava i lable in these references 
' l): lFossati et all (l2007h : [Gil-Hutton & BenavidezI 


(l2003h . There is dependence on telescope focus and sensitiv¬ 
ity to strong wind-induced telescope pointing jitter. We find 
the origin of this continuum polarization variance comes from 
changes in throughput between both polarized beams produced 
by the Savart plate analyzer. Instabilities in the telescope point¬ 
ing, guiding and tracking when combined with intra-order con¬ 
tamination between independent modulated beams induce time- 
dependent continuum changes. In this configuration there is a 
mostly static spectrograph calibration function in addition to the 
induced polarization. As an example, Figure [5] shows repeated 
observations of the unpolarized standard star HR2943 recorded 
witho ut using tip/tilt slit-guiding or charge-shuffling (iGehrelsI 
I 1974 I) . We performed the typical Stokes definition modulation 
scheme for beam-swapping used by most night-time spectropo- 
larimeters. This consists of 6 exposures with the modulators 
tuned such that two modulation states are dedicated to measur¬ 
ing a single Stokes parameter. Figure [5] shows quv in the top 
three panels. The corresponding detected count levels are shown 
in the bottom panel. There are substantial spurious variations in 
quv without any systematic temporal trends. 
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Fig. 6. The computed Stokes u continuum polarization for the unpolar¬ 
ized standard star HD 38393 observed January 10^^ 2013. Spectral aver¬ 
aging by a factor of 80x for high SNR has been applied. At this spectral 
binning there are 50 spectral points per spectral order. Each spectral or¬ 
der is shown in a different color for this wavelength range. The different 
colors nearly overlap but with 0.1% offsets. The intra-order and inter¬ 
order variation is a stable function representing the ~0.5% variation 
across each spectral order. The telescope induced polarization causes 
the curves to change slowly with time and with wavelength. The total 
time elapsed during this data set was 92 minutes (09:59z to 11:31z). The 
statistical noise limits are smaller than 0.05% per pixel at this spectral 
binning. 


When Hi VIS charge-shuffling is enabled, a complete polari- 
metric data set consists of 3 exposures and 6 modulation states 
delivering 12 intensity modulated spectra. When recording two 
modulations states in a single exposure with the new guiding 
system enabled, stable continuum and spectrograph calibration 
measurements are obtained. Figure [7] shows a sequence of 10 
Stokes u polarization measurements obtained over 1.5 hours. 
We observed the unpolarized standard star HD 38393 with 4Hz 
modulation and nominal Hq, Stokes definition liquid crystal volt¬ 


ages (lGehrelsl[l974l) . The standard data reduction scripts were 
used to compute the polarization. The data in Figure [6] have 
been spectrally binned by a factor of 50 to reduce the statisti¬ 
cal (shot) noise limits below 0.05%. A consistent -0.5% varia¬ 
tion remains across individual spectral orders that comes most 
likely from beam throughput variations and intra-order beam 
contamination. This variation we call the spectrograph calibra¬ 
tion function. Each color shows a different spectral order. There 
are discontinuities at orders boundaries from changes in optical 
transmission across a spectral order. However, these discontinu¬ 
ities are now stable and can be calibrated. The computed contin¬ 
uum polarization includes any differential instrumental through¬ 
put variations, since the continuum polarization is calculated us¬ 
ing a simple normalized double-difference: 


Q I a - b c - d 

I 2. G b c d 


( 1 ) 


This normalized double difference method is quite common 
and can easily be compared with other difference and ratio meth- 
ods (iHarrington & Kuhnll2008l:lBagnulo et ^120091 [20121: Patat 
& Romaniello l2006l : lElmore et al .11 19921) . With 4Hz modulation 
and a single-exposure per polarization spectrum, stable varia¬ 
tions with time are seen. Figure [7] shows the residual temporal 
changes for the 10 Stokes u measurements of HD 38393. The 
average u spectrum has been subtracted from each individual u 
measurement. Each color shows a different spectral order and 
there is continuity across order edges. The intra-order variation 
was effectively constant and has been removed. Smooth wave¬ 
length dependent changes of <0.3% per hour over the 1.5 hour 
period are expected as the telescope tracks the star. This contrasts 
with the erratic 5% variation seen in Eigure [5] without charge 
shuffling or guiding. 
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Fig. 7. This Figure shows the residual Stokes u temporal variation for 
the 10 sequential HD 38393 observations from Figureobserved Jan¬ 
uary 10^^ 2013. Each spectral order is a different color. The variation 
with time is a smooth function and all spectral orders overlap for all 
exposures. See text for details. 


We measured the continuum polarization stabilization as a 
function of modulation speed. The dichroic-assisted tip-tilt guid¬ 
ing at 0. IHz was enabled. A single exposure was used to measure 
the continuum polarization spectrum. In order to isolate the opti¬ 
cal stability of the system, the liquid crystal modulators were set 
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to a constant voltage so that all artifacts are caused entirely by 
optical instabilities with minimal impact from real polarization 
changes induced by tracking the star. The target was observed 
for 8 to 12 exposures at a single modulation frequency and the 
root-mean-square (RMS) continuum polarization variation was 
computed. This RMS was measured at modulation frequencies 
of 1/60, 1/15, 1, and 4 Hz with a constant 120 second exposure 
time. The spurious instrumental continuum polarization varia¬ 
tion is 1-4% when doing normal expo sure-time-limited obser¬ 
vations. With modulations at 1-minute scales (l//60Hz) and the 
4 charge-shuffled beams recorded in a single exposure, the con¬ 
tinuum stabilizes to the 0.1% level. When modulation rates are 
increased to the IHz range, continuum polarizations were stable 
to the 0.03% level. 

The optimal modulation speed is a trade off between mini¬ 
mizing spurious effects of optical instabilities and other sources 
of instrumental error. As the modulation speed increases, the sys¬ 
tem spends a larger fraction of the integration time in the process 
of switching states. This reduces the amplitude of the detected 
polarization signal. A key conclusion is that modulations speeds 
faster than IHz show degraded polarization amplitudes. Addi¬ 
tionally, for longer exposures there are more shuffles per expo¬ 
sure. Errors caused by the detector including hot pixels, charge 
traps and charge transfer inefflciency increase. 

2.1. Polarized Standard Stars 

The distribution of polarized standard star reference targets is 
quite uneven and highly clustered at specific locations on the 
sky. It is often not the case that a polarized standard and an 
unpolarized standard are available at the same telescope point¬ 
ing as your science targets in the same season. Spectropolari- 
metric standards are nearly non-existent and interpolation from 
Serkowski-law fit parameters is subject to errors. 

We observed a polarized standard star, HD 30657 with a 
known linear continuum polarization of 3.72% in an R fi lter to 
demonstrate the calibration process (IWhittet et al.lll993) . The 
unpolarized standard star with the closest telescope pointing dur¬ 
ing our allocated observing time was HD 38393. Unfortunately, 
the HD 38393 unpolarized standard star observations were not 
recorded at an identical telescope pointing to our polarized stan¬ 
dard HD 30657. The unpolarized standard star spectra do not 
have any intrinsic polarization and hence record only the tele¬ 
scope induce polarization and spectrograph calibrations. The un¬ 
polarized standard star observations measure the I to QUV terms 
of the Mueller matrix. These terms are nominally stable in time 
and can be subtracted. However, the telescope pointing was not 
identical for the unpolarized and polarized standard star mea¬ 
surements. Thus, the subtracted induced polarization will suffer 
from a small unknown offset. 

The typical change in telescope induced polarization with 
pointing changes was shown in Figure [6] and is fairly small. 
As shown in Figure [7] the u continuum polarization changes 
by roughly 0.2% over roughly one hour of tracking. The aver¬ 
age quv spectra for the HD 38393 unpolarized standard star are 
shown in Figured The spectrum is dominated by u polarization 
and spectrograph intra-order variations. 

We simply subtract the HD 38393 unpolarized standard star 
calibration measurements shown in Figure [8] from observations 
of the polarized standard star HD 30675 taken on the same night 
at a nearby telescope pointing. This allows us to remove the I to 
QUV artifacts within the uncertainty of the different telescope 
pointings. The resulting quv spectra from HD 30675 are thus 
smooth functions of wavelength calibrated for I to QUV effects. 
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Fig. 8. The average quv continuum polarization calculation for the un¬ 
polarized standard star HD 38393 observed January 10^^ 2013. The quv 
spectra are shown over the entire HiVIS wavelength range. The domi¬ 
nant factor is the typical spectrograph calibration in the charge shuffling 
configuration. There were 10 complete measurements recorded at 4Hz 
modulation rates recorded and then averaged. Certain spectral orders 
are shown in color for direct comparison with Figures |7] andSee text 
for details. 


Figure 0 shows these calibrated measurements as the colored 
symbols with red as q, green as u and blue as v. The degree of 
polarization (p) is also shown as the black symbols. The points 
have been spectrally binned by a factor of 50 to decrease the sta¬ 
tistical noise on this V=7.6 star. There is strong cross-talk from 
this linearly polarized standard star as the v spectrum is a nearly 
constant 2% with substantial changes with wavelength. 

The degree of polarization (black symbols) measured for HD 
30675 in Figure 0 closely m atches the 3.72% R-band polariza¬ 
tion of IWhittet et al.l (1 19921) . As the polarized and unpolarized 
star observations were not taken at the same telescope point¬ 
ing, there is uncertainty in the telescope induced polarization. 
There is little variation across spectral orders seen in Figure 0 
The predominant effects seen in continuum polarization is shot 
noise and the telescope induced cross-talk. With these newly- 
stabilized continuum polarization measurements, we can now 
assess our telescope calibrations with much higher accuracy. 

2.2. Standard Star Summary 

We have measured polarized and unpolarized standard stars 
at resolutions of 12,000 and spectral sampling of 100,000 to 
demonstrate calibration of the spectropolarimeter and proof of 
the data analysis algorithms. The calibrations were spectrally av¬ 
eraged to 50 points per spectral order for an effective 1-point 
sampling of R~1250. The calibrations are smooth functions of 
wavelength and show well resolved smooth variation across in¬ 
dividual spectral orders. We have demonstrated stable contin¬ 
uum polarization calibrations at these spectral samplings. Mod¬ 
ulation at speeds of IHz give stable measurements of both po¬ 
larized and unpolarized standard stars below 0.1% precision lev¬ 
els when spectrally binned to high SNRs. This system substan¬ 
tially reduces the types of instrumental effects when comparing 
multiple lines across the spectrum. We find that relatively slow 
modulation rates combined with accurate guiding are sufflcient 
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Fig. 9. The continuum polarization for a polarized standard star HD 
30657 observed on January 10^^ 2013. In Whittet et al. 1992 this star 
is quoted to have an R-band linear polarization of 3.72%. Stokes q is 
shown in red, u in green, v in blue and the total degree of polariza¬ 
tion in black. The colored symbols show the spectrum rebinned to 20 
points per spectral order to highlight noise levels and any order edge 
mismatches. The solid colored lines show the spectra plotted with one 
point per spectral order. The telescope induced polarization and associ¬ 
ated continuum errors have been removed from this measurement using 
an unpolarized standard star observed at a nearby pointing. Observa¬ 
tions were recorded from 11:48 to 12:46. The telescope pointing varied 
from elevation,azimuth of (35,291) to (22,294). Two sets were recorded 
with 300 shuffles per exposure and a IHz modulation frequency. This 
corresponds to 300 seconds per beam, four beams per exposure with a 
total of three exposures. 


to stabilize the optical path of the instrument and remove effects 
from pointing jitter. 

3. Calibrations: Daytime Sky Vaiidation 

Measuring absolute polarization with HiVIS requires a full- 
Stokes polarization calibration process to account for all the tele¬ 
scope and instrument induced cross-talk. In coude path instru¬ 
ments and complex systems, this absolute calibration is often 
the most unstable and difficult. The daytime sky is a very bright, 
highly polarized and easily modeled calibration source. It pro¬ 
vides linear polarization at any telescope pointing, has similar 
surface brightness to bright stellar targets and does not waste 
night time observing. We have developed techniques to observe 
the daytime sky polarization on our low resolution spectrograph, 
LoVIS in order to deri ve polarization calibrati ons for the tele- 
scop e and instrument (iHarrington et al.ll2Qlll). In Harrington 
et al. (l2QI0l) we presented initial calibrations of the HiVIS liq¬ 
uid crystal polarimeter. Here we have successfully applied this 
daytime sky techniques to HiVIS and the liquid crystal demod¬ 
ulation algorithms to our new charge shuffling instrument. Sim¬ 
ple procedures can project these measurements to all telescope 
pointings and calibrate science observations. 

In order to derive the polarization calibrations and the system 
Muller matrix, full-Stokes observations must be obtained of the 
daytime sky polarization spectra over a span of typically a few 
hours. To obtain demodulation matrices, observations of calibra¬ 
tion lamps taken with our polarization calibration injection unit 
are also recorded and processed. This unit consists of a rotating 


Fig. 10. The demodulated daytime sky polarization spectra on May 5th 
2012. The telescope pointing was azimuth 360 and elevation 90. Data 
was taken at two different times (2.62 and 4.86 hours UT). Each mea¬ 
surement has two independent demodulation calibrations applied with 
the two calibration unit demodulation matrices. The curves largely over¬ 
lap showing minimal chromatic effects from the polarization calibration 
unit. Spectral averaging by a factor of 32 was applied to 125 data points 
per spectral order to achieve high SNR. The consistency of the data in¬ 
side each spectral order as well as agreement in regions of spectral order 
overlap is seen in the smooth spectral variation of the quv spectra. 


wire grid polarizer and a rotating achromatic quarter wave plate. 
The measured daytime sky polarization spectra are demodulated 
with our standard procedures using the polarization calibration 
unit optics. Note that for error checking, we record two com¬ 
pletely independent demodulation matrices using our polariza¬ 
tion calibration stage. All polarization calibration unit optics are 
rotated by 180° so that we have two independently calibrated 
demodulated spectra for comparison. 

Examples of demodulated daytime sky polarization spectra 
are shown in Figure [TOl The spectra have been binned to 125 
measurements per spectral order (from an original 4000) in or¬ 
der to highlight systematic errors. The statistical noise is smaller 
than the width of the line used for plotting. One can see in the 
Figure that all curves overlap reasonably well and show only 
small intra-order effects. There are two curves per Stokes pa¬ 
rameter corresponding to the two redundant demodulation ma¬ 
trices. To demonstrate the linear to circular cross-talk induced by 
the telescope, observations at two separate times are shown. The 
sun moved substantially during the 2.2 hours elapsed between 
observations so that the input polarization angle to the telescope 
is known with two values. There is a substantial amount of cir¬ 
cular polarization (blue) measured at short wavelengths caused 
by telescope cross-talk. One can see a substantial variation in 
the measured q and u amplitudes with wavelength indicating the 
telescope induced rotation of the linear polarization orientation. 

In order to derive the telescope calibr ations as a rota¬ 
tion in quv space, the procedures outlined in IHarrington et al.l 
(I 2 OIII) are followed. The degree of polarization is computed as 
^^q^ -vv'^ for each spectrum after demodulation. Each mea¬ 

surement is then scaled to 100% polarization which places the 
measurement on the Poincare sphere. The degree of polariza¬ 
tion expected for simple Rayleigh scattering in the atmosphere 
follows a simple relation whose ra nge of validity has been ex- 
perimentally verified at many sites (ICoulsonll 19801 : ICronin et~^ 
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Fig. 11. The measured degree of polarization for the HiVIS day¬ 
time sky polarization spectra. The telescope was pointed at the zenith 
(elevation=90°) as the telescope azimuth was changed between [90, 
180, 270 and 360] as a sequence in time. The different symbols show 
two different spectral orders: 6700A as diamonds and 8200A as trian¬ 
gles. Different colors represent different telescope azimuths as outlined 
in the legend. As the sun sets, the Zenith degree of polarization is ex¬ 
pected to rise nearly linearly as the solar elevation angle decreases. 


The degree of polarization for a typical set of measurements 
when pointed at the Zenith is shown in Figure [TT] The linearly 
increasing trend in degree of polarization is expected as the scat¬ 
tering angle (sun-zenith-telescope) increases. The slight differ¬ 
ence in polarization between East (090) and North (360) az¬ 
imuths reflects telescope depolarization as well as errors in de¬ 
modulation. The measurements in Figure [TT] are consistent with 
that Rayleigh sky model and a very high (85%) maximum de¬ 
gree of polarization. These values are typical and are similar 
to all-sky polarization meas urements of Haleakala, Mau na-Loa 
observatory, and other sites (iDahlberg et al.ll20TTI.l2009l: Swin- 
dle & Kuhn l2014l) . The high measured degree of polarization 
shows that HiVIS with the liquid crystals and charge-shuffling 
on AEOS has little depolarization, certainly less than 5%. As 
with LoVIS and our other observing modes, observations of a 
wire-grid polarizer covering the calibration sources at the coude 
room entrance show polarization detected above 98% and high 
instrument detection efficiency. 

3.1. Stability and Accuracy of Calibrations 

Our daytime sky calibration technique only requires a minimum 
of two observatimis^to derive a Mupller matrix estimate (Har¬ 
rington et al. l20HT) . Additional observations provide redundant 
information that can be used to test systematic errors and im¬ 
prove accuracy. Independent calibrations can be derived using 
data sub-sets to assess calibration systematic errors. A large 
dataset of daytime sky observations was recorded with HiVIS 
using charge-shuffling between April and June 2012 for verifi¬ 
cation. Calibrations were recorded at cardinal pointings (North, 
East, South and West azimuths) at elevations of 90°, 70° and 55°. 
We collected 4 to 9 measurements at many telescope pointings 


and repeated observations over multiple epochs. For our statisti¬ 
cal analysis, we chose data sets with the telescope pointed at the 
zenith for the cardinal azimuths recorded on April 2V\ 22”^ and 
May 10^^ 2012. For each observation date, an independent sys¬ 
tem M ueller matrix was computed following ([Harrington et al.l 
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Fig. 12. A histogram of the residual rotation angle difference between 
theory and calibrated HiVIS daytime sky polarization spectra. The an¬ 
gular difference between theoretical Rayleigh sky input polarization and 
the calibrated HiVIS measurements on the Poincare sphere is computed. 
There are 19 separate spectral orders spectrally averaged to 50 points 
per spectral order used to create the Histograms. Calibrations were de¬ 
rived using April 21st 2012 data and applied to observations from April 
22nd 2012 giving residual angular error of 1° to 2° between theoretical 
sky models and calibrated observations. Calibrations were also derived 
using both April 21 and 22nd 2012 data which were applied to May 
10th 2012 observations. In this case, the residual angular error is more 
typically in the range of 3° to 6°. 


Since this calibration method effectively applies rotations on 
the Poincare sphere, we found it instructive to show the rota¬ 
tional error between the calibrated observations and the theoreti¬ 
cal daytime sky polarization. First we compute the demodulated, 
fully polarized observations. We then apply our three separate 
system Mueller-matrices to calibrate all the sky measurements. 
This gives us three sets of data calibrated using three indepen¬ 
dent system calibrations. The residual angular distance on the 
Poincare sphere between these calibrated observations and the 
theoretical Rayleigh sky polarization is computed for each wave¬ 
length and telescope pointing. 

Figure [T^ shows the resulting histograms of these residual 
angles. Each histogram curve has been normalized by the total 
number of data points so the relative fraction can be easily seen. 
The solid curve shows a histogram computed using only April 
2V^ and 22"^ data sets and calibrations. When calibrations and 
observations are taken within a 24 hour period, the residual an¬ 
gles are quite small, typically around 1° rotational error in quv 
space. For this distribution, 50% of the points are within 1.3 ° 
and 90% are within 3.1°. When the May 10^^ data is used to 
calibrate April observations and vice versa, the residual error in¬ 
creases. The dashed line in Figure [12] shows calibrations applied 
after twenty days of elapsed time and a reinstallation of the po- 
larimeter. This histogram shows a peak around 4° rotational er¬ 
ror. This distribution has 50% of the points within 3.9° rotational 
error however 90% of the points are still within 6.5°. 
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There are several potential causes for calibration error. First, 
the Rayleigh sky theory can show angular variations between 
actual observations at high precision. Most all-sky imaging po- 
larimetry studies have focused on the degree of polarization but 
the angle of polarization is not monitored as precisely (Voss & 
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Recent experiments with a new all-sky imaging polarimeter on 
Haleakala adjacent to AEQS have inve stigated this uncertainty 
([Swindle & Kuhnll20l4 [Swindle![2014]) . They find that the an¬ 
gle of polarization can vary by up to 3 degrees in regions of the 
daytime sky with degree of polarization lower than about 15%. 
Provided we only use daytime sky measurements in regions with 
measured degree of polarizations above 15%, this angular varia¬ 
tion is smaller than 1° in input qu orientation. Thus, we only use 
daytime sky measurements with high degree of polarization and 
minimize this source of uncertainty. 

Thermal stability of our liquid crystal retarders is another er¬ 
ror source. Though we are in a coude room, small changes in 
temperature during the time between modulation matrix mea¬ 
surement and sky observation can cause errors. Self-heating 
from the control signal is also known to cause retardation 
changes though this effect is small in our devices. A n exam¬ 
ple can be seen in Figure 5 of [Harrington et al.[ ([201Qh . A small 
temperature change between April and May could easily cause 
an instability. Replacement of our liquid crystals with the active 
temperature stabilization version and better thermal control of 
the instrument could possibly improve our calibration stability. 



Fig. 13. This Figure shows the six Mueller matrix estimates derived 
via our least-squares method using all data from our 2013 to 2014 cal¬ 
ibration efforts. Each colored point represents an altitude-azimuth tele¬ 
scope pointing where we have quality data for estimating the telescope 
Mueller matrix. The QQ and QU terms are on top. The UQ and UU 
terms are in the middle. The QV and UV terms are on the bottom. Note 
the bottom row represents linear to circular polarization cross-talk and 
the amplitudes are at 100% for many pointings. 


3.2. Altitude-Azimuth Dependence 

Telescope polarization calibration must be applied to observa¬ 
tions at all telescope pointings. Creating a telescope Mueller 
matrix at every altitude and azimuth requires sufficient calibra¬ 
tion data to accurately characterize the behavior of the telescope 
Mueller matrix. An observing campaign was performed during 
winter of 2013 and summer of 2014 to make a densely covered 
set of calibrations with Hi VIS and charge shuffling. Daytime sky 
polarization observations were obtained in October, November 
and December of 2013 as well as in May of 2014. We recorded 
many redundant exposures at azimuths of [030, 060, 090, 120, 
150, 180, 210, 270, 330, 360]. The corresponding elevations 
were: [10, 20, 25, 35, 50, 60, 75, 89]. Recent upgrades to the 
AEOS telescope now allow observations all day. Previously, our 
telescope pointing elevation limit was above 55° and the sun had 
to be below elevations of 45°. 

The procedure outlined in [Harrington et al.[ ([2011[) allows one 
to estimate Mueller matrix elements from Rayleigh sky obser¬ 
vations. We scale all HiVIS daytime sky observations to 100% 
degree of polarization before we calculate the cross-talk terms 
of the Mueller matrix. Measurements are put on to the Poincare 
sphere before deriving the best-fit rotation matrices as this re¬ 
sults in constructing a physically realistic Mueller matrix. This 
procedure models telescope as a simple rotation incoming quv 
polarization to some other quv basis set. As an example of the 
altitude-azimuth behavior, we show the Mueller matrix estimates 
for spectral order 3 in Figure [131 

For AEOS, the functional dependence with altitude and az¬ 
imuth follows simple trigonometric functions and can be fit with 
a few independent parameters. The dominant polarization cross¬ 
talk is caused by rotation of the fold mirrors along the altitude 
and azimuthal axes. With simple functional coefficients, we can 
re-compute the azimuthal dependence at an arbitrary telescope 
pointing. Figure [T7] shows the result of this interpolation. The 
cross-talk elements of the telescope Mueller matrix are shown 
for all altitudes and azimuths. This process allows individual 
Mueller matrix elements to be obtained at the pointing of the 
astronomical target. After interpolation, the user can derive the 
best-fit rot ation matrix using the stp dard Euler-angle procedure 
outlined in [Harrington et al.[ ([201 ih . This ensures physicality of 
the telescope Mueller matrix and reduces errors in Mueller ma¬ 
trix term estimation. By representing the telescope Mueller ma¬ 
trix as a rotation matrix, we ensure physicality of the transfor- 
mation and we do not amplify certain kinds o f systematic noise 
([Kostinski et al.[[l99^:[Takakura & Stoll2009[: Givens & Kostin- 
ski M993h . 

To quantify the errors, we derive the cumulative distribution 
function for the difference between the six independent Mueller 
matrix estimates and the corresponding best-fit rotation matri¬ 
ces after applying the altitude-azimuth interpolation. Figure [14] 
shows each Mueller matrix estimate in a different color. The 
trigonometric function fitting results in an uncertainty of less 
than 0.1 between the estimates and the rotation matrix fits. 

One of the major sources of calibration uncertainty in the 
present approach is the stability of the liquid crystals in tem¬ 
perature and time. The individual measurements can calibrated 
and differenced from the theoretical daytime sky polarization to 
demonstrate the demodulation and cross-talk calibration errors. 
Figure [15] shows the difference between the theoretical Rayleigh 
sky polarization spectra and demodulated calibrated HiVIS mea¬ 
surements. In this Figure, April 2 2012 data was used to de¬ 

rive Euler angles. This data was then polarization calibrated and 
subtracted from the theoretical Rayleigh sky polarization. The 
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Fig. 15. This Figure shows the difference between the demodulated and calibrated Stokes parameters and the predicted Rayleigh sky polarization 
spectra for April 2012. Each panel shows observations calibrated at a different telescope azimuth. There are 4 to 5 independent observations 
at each telescope pointing. Different colors represent residual Stokes parameters for individual observations computed as data minus model. The 
residuals are always below 0.1 across the wavelength range. Note - only a single calibration point per spectral order has been applied and the data 
is shown binned to 50 spectral samples per order. Note that these residual Stokes parameters give rise to the few degree rotational error between a 
measured Stokes vector and a model Stokes vector projected on the Poincar’e sphere as shown above in Figure[T2] 



Mueller Matrix Element Fit Error 


Fig. 14. The cumulative error distribution functions when fitting rota¬ 
tion matrix functions of azimuth and elevation to the individual mueller 
matrix elements. The distribution is based on the difference between the 
best-fit rotation matrix and each individual least-squares based Mueller 
matrix estimate at every measured azimuth and elevation. Each color 
corresponds to one of the six Mueller matrix estimates computed us¬ 
ing our daytime sky method. The 68% confidence interval falls between 
errors of 0.04 for some elements and 0.10 for other elements. 


residual errors are typically less than 0.05 in any individual 
Stokes parameter. These residual errors were used to derive a 
rotational measure of absolute calibration accuracy and are the 
basis for the histograms of the angle between the measured cali¬ 
brated Stokes vector and the theoretical daytime sky Stokes vec¬ 
tor shown above in Figure [121 

We find that the likely cause of the ripples seen at small 
wavelength scales within individual spectral orders is caused by 
slight changes in liqu id crystal^properties. In Figure 3 of Har¬ 
rington et al. (120101) 7 we showed a transmission function using 


a very small beam footprint using a Varian spectrophotometer. 
This Varian is effectively a lab-calibrated commercial scanning 
monochrometer. This instrument also showed small-scale ripples 
in the overall transmission function that were highly dependent 
on placem ent of the liquid crystal. Similar ripples were seen in 
Figure 5 of lHarrington et al.l ( |2_Q1Q|) when the temperature of the 
device was changed. Liquid crystal stability will likely be a lim¬ 
itation of our system. We note that the calibrations applied here 
only use one telescope Mueller matrix element per spectral or¬ 
der. Potentially performing calibrations with more data at higher 
spectral sampling could improve the calibration procedures. 

The wavelength dependence of the telescope Mueller matrix 
is very slow. Figure [16] shows the six Mueller matrix estimates 
after being ht to simple trigonometric functions of altitude and 
azimuth. The top panel shows the 6225A spectral order while 
the bottom panel shows the 8601A spectral order. The linear po¬ 
larization terms are dominated by the transformation between an 
altitude-azimuth based frame and a Zenith frame. The circular 
polarization terms are quite strong at 6225A but tend towards 
zero at 8601A. Given this gradual wavelength dependence, we 
hnd that one calibration per spectral order is sufficient for this 
demonstration. 


3.3. Example Calibration of e Aurigae 


Epsilon Aurigae is a very bright tarj 
variable Hq, polarization signature ( 

?et with a known strong and 

E 

[arrineton & KuhnI 2009al: 

iKemp et al.lll986 

;ICole & Stencel]l20 

lll:lKemDetal.lll985tHall 

& Henson I 2 OIOI; 

Henson et al.ll2012l; 

Burdette & Hensonll2012l; 


IStencelll2013l) . We have a long and detailed set of e Aur observa 
tions taken with ESPaDOnS at CFHT for comparison. We have 
never detected any circular polarization in Hq, in the 2006 to 2012 
time frame. As an example of this stars Hq, line prohle, we com¬ 
piled 42 independent measurements in to blocks where the target 
shows substantially different polarization and intensity line pro- 
hles. Figure [18] shows the observations broken in to 9 main ob¬ 
serving blocks covering summer and winter observing seasons 
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Fig. 16. This Figure shows the six Mueller matrix estimates trigonom¬ 
etry function fitting for all altitudes and azimuths. The QQ and QU 
terms are on top. The UQ and UU terms are in the middle. The QV 
and UV terms are on the bottom. The top graphic shows spectral or¬ 
der at 6225A. The bottom graphic shows spectral order at 8601A with 
much reduced circular polarization cross-talk (low QV and UV values). 
The sinusoidal variation in the linear polarization terms represent the 
transformation between the altitude-azimuth frame and a Zenith based 
polarization frame. 


over the 6 years. We followed our recipe for post-processing of 
the Libre -Esprit reduced ESPaDQ Ns spectra outlined in the ap¬ 
pendix of lHarrington et al.l (l2009l) . Essentially our process com¬ 
bines the reduced spectra from individual spectral orders and av¬ 
erages pixels to a more uniform spectral sampling. In the case 
of Eigure [18] the effective spectral sampling is at about 35pm 
per pixel for a 1-point sampling of R~ 170,000 data. Eor the pur¬ 
poses of this calibration we note that the polarized Hq, profiles 
have amplitudes of 0.5% up to 1% while Stokes v signatures are 
undetectable to within the signal-to-noise thresholds. The effec¬ 
tive signal to noise ratios in both the ESPaDOnS data sets is over 


1,000 per pixel and the relevant spectral features are very broad. 
Our Hi VIS charge-shuffled spectra are sampled at ~44pm per 
pixel at Hq, and we highly over-sample the delivered optical reso¬ 
lution of 12,000 to achieve high signal-to-noise ratios efflciently 
on bright targets. 

During October 23, 2013, a large set of observations were 
recorded on the star e Aurigae while performing our daytime 
sky calibration mapping procedures. We collected 30 individual 
polarization measurements over the course of one night. Over the 
observations, azimuths ranged from 055° through transit to 330°. 
Elevations were between 34° and 67°. Each full-Stokes e Aur 
observation set was reduced and demodulated using the standard 
observing sequence with our polarization calibration unit. The 
signal-to-noise ratio for each individual quv measurement was 
between 650 and 900 for all 30 independent quv measurements. 



6555 6560 6565 6570 
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Fig. 18. This Figure shows the 6 Aurigae Hq, line profiles observed 
with ESPaDOnS between 2006 and 2012. This star has shown strong 
variable qu profiles but Stokes v was always zero within the statistical 
noise. See text for details. 


There were 4 separate clusters of data recorded at telescope 
pointings of (051,35), (045,52), (011,67) and (333,64). There 
are 5, 10, 5 and 10 consecutive exposures in each of the associ¬ 
ated observing blocks. Eigure [19] shows the quv spectra averaged 
during these 4 observing blocks. The ESPaDOnS spectra show 
that the Hq, signature is always contained entirely in linear po¬ 
larization qu. However, since the linear to circular cross-talk is 
100% at some telescope pointings for Hi VIS, our uncalibrated 
data shows this Hq, signature largely in Stokes v. Some changes 
can be seen in the details of the quv profiles, showing the chang¬ 
ing cross-talk as the telescope tracks the star. 

Using the telescope Mueller matrix element estimates, the 
best-fit rotation-matrix was computed for the telescope pointing 
at each e Auri gae data set f ollowing our procedures (Harring¬ 
ton et al. 120111) . The telescope Mueller matrix maps of Eigure 
[TJcan be interpolated to the telescope pointing of every individ¬ 
ual observation. Eigure [20] shows the telescope Mueller matrix 
estimates for every individual observation of e Aurigae. The 4 
observing blocks can be seen as clusters of points with similar 
Mueller matrix elements. Eor example, the far right side of Eig¬ 
ure [20] shows a cluster of 5 observations while e Aurigae was 
rising at (azimuth, elevation) of (051,35). The telescope point¬ 
ing change during this time was largely elevation. On the far left 
side of Eigure [20] a sequence of 10 observations is seen near an 
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Fig. 17. An altitude-azimuth map of the system Mueller matrix cross-talk elements for the 6560A spectral order. The best fit rotation matrix was 
fit to the Mueller matrix estimates elements following our standard procedure. The Mueller matrix element az-el maps are all linearly scaled from 
black at -1 to white at -i-l. Each element is plotted from the horizon to the Zenith for the full 360° azimuth range. 


(azimuth,elevation) of (333,64) while the predominant telescope 
pointing change was in azimuth. 

These telescope Mueller matrix calibrations were applied to 
each 6 Aurigae data set. Figure ^T\ shows the average demodu¬ 
lated and calibrated data sets for the 4 major observing blocks. 
The top three panels show Stokes q, u and v respectively. The 
bottom panel shows the continuum-normalized intensity profile 
for Hq,. The uncalibrated measurements showed the observed 
signature was dominantly Stokes v with substantial q and u. Af¬ 
ter applying the telescope calibrations, the Hq, line profile signa¬ 
ture is almost entirely contained in Stokes q and u. 

There are several sources of error in these demonstration 
calibrations that can easily be improved. First, we have only 
used one calibration per spectral order which gives some mild 
wavelength dependent residual errors across a spectral order. 
Secondly, the daytime sky observations were recorded over a 6 
month time frame. This leads to several degrees rotational er¬ 
rors as seen in the histogram of Figure [121 Note that a Stokes 
parameter error of 0.05 corresponds to a 3° misalignment while 
a 0.2 error corresponds to an 11° misalignment. Since we only 
apply one calibration per spectral order, for an effective reduc¬ 
tion in wavelength sampling by a factor of 4000, we get resid¬ 
ual errors wavelength from the liquid crystal ripples and general 
wavelength dependent behavior as seen in Figure [TSl The third 


error comes from interpolating the best-fit rotation matrix to the 
telescope pointings of every individual observation. The cumu¬ 
lative distribution functions of Figure [14] show that the Mueller 
matrix element fitting process gives Icr errors at element values 
of roughly 0.04 to 0.10 for a 68% error level. This corresponds 
to additional rotational error levels of 2° to 6°. 

The calibrated quv behavior for e Aurigae is shown in Fig¬ 
ure [^ To show the quv behavior at high signal to shot-noise 
ratios, we combined all individual calibrated spectra and then 
averaged 4 pixels in the spectral direction. The signal to noise 
ratio of these average spectra is between 5700 and 7700 as esti¬ 
mated from the standard deviation of the noise of the continuum 
adjacent to the Hq, line. After the 4-pixel averaging, we sample 
our R~ 12,000 spectra at an effective sampling of 34,000 for 3 
points per resolution element. The residual Stokes v signature 
in Figure [^ is expected is an amplitude of 0.14% while the qu 
amplitude is at 0.50%. This corresponds to a misalignment of 
15.6°. As shown in Figure [T^ we expect leakage of this magni¬ 
tude considering all the instrument changes in the long time span 
between calibrations and observations. 

4. Conclusions 

We have demonstrated here general techniques that overcome 
both calibration-type uncertainties (I) and measurement stabil- 
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Fig. 19. This Figure shows the average quv data sets in 4 main telescope 
pointing blocks for e Aurigae observed with HiVIS on October 23^^ 
2013. No telescope Mueller matrix calibration has been applied. The 
input signal is qu as shown by ESPaDOnS observations but the detected 
HiVIS signature is much larger in v due to the telescope cross-talk. The 
4 independent quv spectra show some spectral variation due to changes 
in the telescope cross-talk as the optics change with time. See text for 
details. 
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Fig. 21. This Figure shows the average of calibrated quv data sets or 
Epsilon Aurigae observed with HiVIS on October 23^^ 2013. The 4 
main observing blocks at similar telescope pointings were averaged to 
achieve high SNR. The signal is largely qu as expected from ESPaDOnS 
observations. The Stokes v measurements are nearly zero. Some resid¬ 
ual Stokes V is expected when using calibrations derived over a 6 month 
time frame. See text for details. 
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Fig. 20. The estimated system Mueller matrix elements changing with 
time for the 30 6 Aurigae data sets collected with HiVIS on the 23rd 
of October, 2013. This shows the time dependence of the system cross¬ 
talk elements of the Mueller Matrix varying quite substantially with 
telescope pointing and time. See text for details. 
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ity issues (II). SNRs above 0.01% were achieved at full spectral 
resolution using calibration sources. On-sky calibrated observa¬ 
tions of 6 Aurigae achieved SNRs of 5000 to 7600 at R=34,000 
sampling. Using our new daytime sky calibration technique, we 
have independently calibrated the polarization response of the 
instrument and telescope as functions of altitude and azimuth 
in a system with 100% cross talk. In this altitude-azimuth tele¬ 
scope the polarization cross-talk calibration has a simple func¬ 
tional form on-sky and can be done easily on a routine basis to 
assure calibration stability. The cross-talk calibration was tested 
on a stellar target e Aurigae with good results. Several sources of 



Fig. 22. This Figure shows the quv behavior of calibrated quv data sets 
or 6 Aurigae observed with HiVIS on October 23^^ 2013. The 4 main 
observing blocks at similar telescope pointings were combined after cal¬ 
ibration to achieve high SNR of 5,000 to 7,600 in quv. The residual er¬ 
rors can be seen as a cluster of points around 0 with cr < 0.05%. The 
signal is largely qu with a linear polarization amplitude of 0.50% as 
expected from ESPaDOnS observations. The calibrated Stokes v ampli¬ 
tude is roughly 0.14%. This residual Stokes v corresponds to a 15° ro¬ 
tational calibration error between linear and circular polarization. This 
error is expected when combining all the uncertainties in our calibration 
demonstration derived over a 6 month time frame. See text for details. 


calibration algorithm uncertainty have been identified and future 
improvements are expected. 

The quality of our cross-talk calibration method was shown 
to be better than 4% in an individual Stokes parameter if cali¬ 
brations are taken within a day of the observations at the same 
telescope pointing. This corresponds to a 2° Poincare-sphere 
rotational uncertainty of the polarization orientation. Calibra¬ 
tions taken even a month after observations still provide 5° rota¬ 
tional accuracy which corresponds to about 10% of an individual 
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Stokes parameter. Our assumptions of a weakly depolarizing our 
telescope and instrument are reasonable given the degree of po¬ 
larization measured for polarized standard stars as well as for 
the high measured daytime sky degree of polarization maximum 
i^max) • 

We were able to measure both polarized standard stars and 
unpolarized standard stars with smooth reproducible continuum 
polarization at the 0.1% level even in the presence of a >2% 
inter-order artifact. We applied firmware and software upgrades 
to accomplish charge shuffling with a CCD that was not origi¬ 
nally designed for this purpose. Many of our instrument stability 
issues (type II) were overcome while improving our measure¬ 
ment efflciency. This new combination of charge shuffling, liquid 
crystal modulation and daytime sky calibration can deliver cali¬ 
brated, high precision observations for use in a variety of stellar 
spectropolarimetric applications when high signal-to-noise data 
are obtained. 

For other night-time instruments there are some relevant 
lessons learned from this study. Often, the major system limi¬ 
tations in high signal-to-noise applications is the ability to col¬ 
lect photo-electrons rapidly and suppress instrumental errors. We 
found that bidirectional CCD charge shuffling yielded signifi¬ 
cant calibration and efficiency gains. Multiple beams and storage 
buffers can be used without readout, reducing overheads. Mod¬ 
ulation can be performed much faster than the exposure time. 
Operational flexibility is easily controlled via software. 

Optical stability is also an important improvement. The ben¬ 
efits of improved guiding, reduced slit height, IHz modulation 
and having complete polarimetric information within a single ex¬ 
posure allowed sensitive continuum polarization measurements 
with the HiVIS spectropolarimeter. Fast (kHz and higher) mod¬ 
ulation are not essential to substantially remove many sources of 
instrument instability and pointing jitter. These optical improve¬ 
ments will also remove sources of calibration uncertainty. 

The daytime sky calibration was readily adapted to this in¬ 
strument. The telescope Mueller matrix elements de-rotated ob¬ 
servations to a Poincare rotational accuracy of a few degrees with 
no night observing time lost. For AEOS, the altitude-azimuth 
dependence of the telescope Mueller matrix has a relatively 
simple functional form, which would be even simpler in equa¬ 
torial or non-coude optical systems. Calibrations at relatively 
few telescope pointings are sufficient to capture most of the ob¬ 
served Mueller matrix element variation. As night-time and solar 
telescopes adopt polarimetric instruments in non-optimal opti¬ 
cal configurations we believe the calibration strategies described 
here will find widespread applications. 
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Appendix A: Scattered Light as poiarization 
contamination 

The optical performance and stability of a spectrograph impacts 
our ability to calibrate systematic errors and to measure sta¬ 
ble continuum polarization. Optical imperfections cause scat¬ 
tered light, wavelength instabilities and spectropolarimetric is¬ 
sues with amplitudes above 0.1%. These are also particularly 
difficult for intra- and inter-order stability which effect abso¬ 
lute polarization calibration. In order to estimate these effects, 
we performed a thorough optical stability analysis and detector 
characterization. 

Flat-field calibration images were computed and averaged to 
very high SNR to measure and remove the scattered light con- 
taminati on. We chose to use the LoVIS configuration I Harring- 
ton et al. ]2^d) with single spectral order imaged on the detector. 
Most of the CCD focal plane is not directly illuminated in this 
configuration allowing us to determine the wide-angle scattered 
light properties without inter-order contamination from adjacent 
spectral orders. Figure lAJl shows the spatial profile derived from 
a 100-exposure average of fiat field measurements scaled to the 
peak flux. The dark and bias values have been removed from 
this stacked image with another 100-image combination of dark 
frames. This ensures high precision characterization of the scat¬ 
tered light in the spectrograph as the dark and bias are indepen¬ 
dently removed without choosing a specific region for dark sub¬ 
traction. The red vertical line shows the edge of the spectral order 
as defined by the data reduction pipeline in the typical HiVIS 
charge shuffling mode. We use a 13 pixel half-width from the 
central spatial pixel for a total width of 27 pixels used in the 
spectral extraction. The flux levels at this order edge are typi¬ 
cally below 1% of the peak brightness. The blue lines show the 
region of the image typically used for background subtraction in 
the HiVIS data reduction pipeline. This region is usually another 
factor of 5 or more below the peak flux levels. Similar results 
for scattered light are seen in the HiVIS configuration with both 
the achromatic rotating retarders and the liquid crystal variable 
retarders. 

The efficient subtraction of scattered light is important for 
absolute polarization calibration. Any background intensity not 
effectively removed from all four polarized beams for each spec¬ 
tral order will depolarize and add other continuum polarization 
effects. There are 17 spectral orders extracted with HiVIS using 
the blue cross-disperser. There are two CCDs with two readout 
amplifiers per CCD. Each spectral order has pixels outside the 
illuminated region that are used for extracting an effective back¬ 
ground. These backgrounds show scattered light at the few hun¬ 
dred count level. 

Additional information about optical contamination in the 
presence of strong and spectrally narrow features is easily de¬ 
rived from arc lamp spectra. HiVIS has a Thorium-Argon (ThAr) 
lamp mounted on a calibration bench at the entrance to the coude 
room. An example of a high-precision monochromatic slit im¬ 
age derived from ThAr exposures is shown in Figure IA.21 The 
red cross disperser was used in HiVIS with the 1.5" slit. No 
charge shuffling was performed to illustrate the optical perfor¬ 
mance at high scattering angles. HiVIS has a spectral sampling 
of 4.5pm to 5.5pm per pixel corresponding to a 1-pixel sampling 
of R~ 150,000. The spectrum is highly over sampled. The Figure 
is log-scaled and one can easily see there are some asymmetries 
in the background flux in the spatial direction at the 0.1% level. 
This is quite small and stable with minimal influence on the de¬ 
livered polarimetric calibration accuracy. 
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Fig. A.l. This Figure shows the spatial profile delivered for the spec¬ 
trograph in the LoVIS configuration with the fiat-field lamp (April 26th 
2011). There were no polarimetric components present in the beam. 
The spatial profile is shown on a log scale. The red vertical line shows 
the edge of the extracted spectrum typically set in our data reduction 
pipeline. The blue lines show the region typically used to compute the 
background (dark, bias & scattered light) extraction. 



80 100 120 140 160 180 200 
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Fig. A.2. This Figure shows a log-scaled extraction of Th-Ar arc line 
profiles computed using a stack of 40 individual exposures from July 
25th 2012. A clear asymmetry is seen in the wings of the monochro¬ 
matic slit-image at roughly 0.01% of the peak intensity. 


This instrument has been used to measure spectropolarimet- 
ric profiles in highly distorted emission lines. Differential bleed¬ 
ing of intensity from different spectral locations can contaminate 
observations. Spectral lines from P-Cygni and Herbig Ae/Be 
stars can have intensity contrasts of >30 over very narrow spec¬ 
tral regions. The effective spectral resolution is typically derived 
by measuring the full-width-half-maximum of arc lines. How¬ 
ever, at small signal levels a more complete characterization is 
necessary. To estimate potential leakage from different spectral 
orders, we computed the spectral profile for arc lamp lines using 
stacks of 40 images. These profiles show that we deliver >99% 


concentration of monochromatic light within 30 spectral pixels. 
The differential contamination for Hi VIS is quite small. 

In addition to spectral bleeding, any Savart plate can poten¬ 
tially cause degradation. A Savart plate is manufactured so that 
the extraordinary and ordinary beams interchange between the 
two crystals and the two exit beams correspond to ^ - 6> and o-e 
pathways. However, imperfect construction can cause leakage in 
to the ee and oo beams We had a special alignment procedure 
done at the manufacturer with a <10“^ leakage specification. 
We also have done tests of our upgraded Savart plate and we 
find this beam leakage undetectable. This specification is critical 
to HiVIS as leakage would have directed two additional spec¬ 
tral beams in between polarized orders decreasing contrast and 
adding spectral leakage at ~54 pixel separations. 

The spatial width of the extracted spectra combined with the 
slit dekker width and any optical imperfections will cause intra¬ 
order contamination between modulated beams. This will cre¬ 
ate a systematic continuum polarization error when doing charge 
shuffling. The light that bleeds from one modulation state with 
polarization (a - bla + b) in to the other modulation state with 
polarization {c - die + d) asymmetrically will cause continuum 
variations. As an illustration of this effect, Fisure I a 31 shows the 
four beams for the 10^^ spectral order for HiVIS on May 25^^ 
2010. In order to exaggerate the typical continuum effects, the 
slit dekker was set wider than nominal. The modulators (liquid 
crystals and rotating retarders) were removed to ensure that scat¬ 
tered light was the dominant error. 



Spotioi Pixel 

Fig. A.3. This Figure shows the computed spatial profile for the spec¬ 
tral order 10 with a wide dekker setting. This wide dekker exaggerates 
the intra-order contamination. There were no modulators in the beam. 
Charge shuffling of 27 pixels was done with polarization calibration 
unit polarizer in place. Beams a and b represent one set of orthogonally 
polarized spectra. Beams c and d represent another set of orthogonally 
polarized spectra. The asymmetry in intra-order contamination is eas¬ 
ily seen where intense d and b beams strongly influence c but not a. 
In double-differenced polarization calculations, this leads to spurious 
continuum polarization effects. 


One can see in Figure IA.3I that the asymmetric bleeding 
of one beam in to the other will cause different computed in¬ 
tensities even though no modulation or optical movement was 
present. The incomplete separation of beams lead to a system¬ 
atic decrease in the detected degree of polarization for each of 
the HiVIS orders as well as an increase in the systematic er- 
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ror across each individual order. The b and d beams had similar 
intensity and similar bleeding so these spectra were influenced 
little. The c spectra was much more substantially influenced be¬ 
cause of it’s low intensity (300 counts) and the proximity to the 
bright h spectrum. The uneven disruption of the a and c spec¬ 
tra caused a much more substantial error in the (c - J)/(c -r d) 
polarization measurement. 

In typical double difference calculations, the two methods do 
not give the same continuum polarization. The two methods in 
this nomenclature are: 


^ ab-cd 


Q 

I 


1 a-b c-d 

2 a + b c + d 


and 


(A.l) 


^ac-bd 


Q 

I 


1 a-c b-d 

2 a + c b d 


(A.2) 


Without asymmetries in the optics, detector and modulation, 
these two e quations should give identical results. In Harrington 
& Kuhn (I2Q08I) this is exactly what we And for HiVIS with only 
two well separated beams. The qab calculation does a normalized 
ratio with spectra recorded simultaneously in time. The qac cal¬ 
culation is a normalized ratio with spectra recorded on the same 
pixel for gain independence but with temporal changes between 
beams possible. If these four beams are not well separated and 
stable in time and completely unpolarized at the analyzer, con¬ 
tinuum polarization errors will be present. As a demonstration. 
Figure |A3] shows the computed polarization for this May 2010 
experiment with the flat field calibration lamps. Since the mod¬ 
ulators were removed and no optical motion occurred between 
modulation states, the derived Stokes q polarization should be 
nearly 100% and identical between the ab and cd calculations. 
Figure IA.4I shows that the ab spectra have a very different con¬ 
tamination than the cd spectra across individual spectral orders 
at the 0.5% level. 

These asymmetries combined with typical optical instabili¬ 
ties (pointing, guiding, seeing, etc) conspire to add systematic 
continuum variations across spectral orders that are difficult to 
calibrate. This is a major limitation when stabilizing the contin¬ 
uum polarization of a slit spectropolarimeter. 
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